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1^ I We study by kinetic Monte Carlo simulations the dynamic behavior of a Ziff-Gulari-Barshad (ZGB) model for the 
. catalytic oxidation of CO on a surface. It is well known that the ZGB model presents a continuous transition between 
an oxygen poisoned state and a reactive state that it is not observed in nature. Based on some experimental results 
, that indicate that the oxygen atoms move away from each other upon dissociation at the surface, we modify the 

■ standard ZGB model by changing the entrance mechanism of the oxygen molecule. We study three different ways 
\ in which the oxygen atoms can be adsorbed at the surface such that the nonphysical continuous phase transition 

■ disappears. We calculate the phase diagram for the three cases and study the effects of including a CO desorption 
' mechanism. 

(-| ; I. INTRODUCTION 

D . Catalytic reactions on surfaces have attracted a great deal of interest. Besides their enormous technological and 
industrial applications, they present a rich and complex behavior that makes them an ideal laboratory to study 
^ . non-equilibrium phenomena [lH3[. 
c5 ■ The Ziff-Gulari and Barshad (ZGB) model 0, is a simple but fruitful system to describe the catalytic oxidation of 
carbon monoxide (CO) on a surface. In this model the catalytic surface is represented by a square lattice whose sites 
"T ■ can be empty or filled with oxygen (O) or CO. When CO and O are located on nearest-neighbors sites on the lattice 
' I they react and liberate CO2, leaving behind two empty sites on the surface, i.e., a Langmuir-Hinshelwood mechanism. 
The kinetic of the reaction is described by a single parameter, y, the probability that the next molecule arriving to the 
surface is a CO. The model presents two phase transitions, a continuous one at j/i between a passivated phase where 
the surface is completely filled with oxygen and a reactive phase, and a discontinuous one at y2 between the reactive 
phase and a passivated phase with the surface filled with CO. However, experimental results indicate that above a 
critical temperature, the discontinuous transition disappears and the CO2 production decreases continuously 

■ This behavior can be reproduced by modifying the ZGB model to include a CO desorption rate, k, that mimics the 
\ effect of the temperature (it has been shown experimentally that the desorption rate of the oxygens is negligible 
. [lo|). Experiments also indicate that in real systems there is not a passivated phase where the surface is poisoned 

with O, the reaction rate increases as soon as y departs from zero [lO: Ell- More realistic treatments of the catalytic 
, oxidation are oriented toward modifying the ZGB model to eliminate the nonphysical transition at yi . Several authors 
have shown that this transition can be eliminated by adding to the Langmuir-Hinshelwood an Ele v-R ideal step that 
allows a reaction between CO molecules in the gas phase and adsorbed O atoms on the surface [13 - ll4| . However, 
there are no experimental data indicating that such a reaction between free CO and adsorbed oxygen occurs. Inspired 
by experiments based on scanning tunneling microscopy that point toward a "hot atom" mechanism according to 
which the two oxygen atoms are propelled apart upon dissociation [Tsf , we show that the nonphysical transition at 
yi can be eliminated by making simple changes to the entrance mechanism of the oxygen atoms in the ZGB model. 
Similar mechanisms have been suggested by other authors [l6l - [T8| . In the ZGB model upon adsorption the oxygen 
atoms are located one lattice unit apart. In this work we compare three different entrance mechanisms for the oxygen 
atoms. In the first case, upon dissociation the oxygen atoms are separated by lattice units, in the second they 
are located at opposite sites of an empty site, separated by 2 lattice units, and in the third case their locations are 
the same as in the second but the center site does not need to be empty. All these mechanisms have the advantage 
of incorporating the hot atom effect without changing the nature of the model. In this work we perform Monte Carlo 
simulations to calculate the phase diagrams of the different models. Then, in order to include temperature effects, we 
add a desorption probability for the CO molecules and determine how the phase diagrams are affected. 
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II. THE MODELS AND THE NUMERICAL SIMULATION 



The models are simulated on a square lattice of linear size L that represents the catalytic surface. A Monte Carlo 
simulation generates a sequence of trials: A CO or O2 molecule is selected with probability y and 1 — y respectively 
These probabilities are proportional to the partial pressures of the gases. The algorithm works in the following 
way. A site i is selected at random. If i is empty and CO adsorption is selected (probability y), a CO molecule is 
adsorbed. If one of the nearest-neighbors (nn) of i is filled with an O, there is a reaction between the CO and the O: 
a CO2 molecule is liberated leaving two empty nn sites at the surface. The O2 adsorption requires two empty sites 
on the lattice to accommodate the two dissociated atoms. There are several ways to choose these sites. The main 
purpose of this work is to compare these different adsorption mechanisms. In the standard ZGB model, if the site i is 
empty and O2 adsorption is performed (probability 1 — y), a nearest neighbor of i is selected at random, if it is empty 
the adsorption takes place: one O is adsorbed at i and the second at the chosen nn site. Another option is to adsorb 
the two oxygens at next-nearest neighbors (nnn) sites: the adsorption proceeds only if one of the nnn of the empty 
site i is also empty, this case is going to be labeled the ZGB-n model. A different mechanism that we are going to 
label the Hal model, consists in performing the adsorption only if the selected site i is empty and one of the vertical 
or horizontal nn pairs of the site i are empty. If this is the case an O is located in each site of the pair (leaving the 
site i empty). In a slightly different version, that will be labeled Ha2, the site i does not need to be empty for the 
adsorption to proceed in one of its nn pairs. For all the models, after an O2 adsorption has proceed, all the nn of 
the adsorbed oxygens are examined. If any of them is filled with an CO a reaction proceeds, and a CO2 molecule is 
liberated. The different adsorption mechanisms for O2 are sketched in Fig. [1] and Fig. [2j By locating the adsorbed 
oxygens further apart than in the original ZGB model, in a certain way all the new mechanisms take into account the 
hot-atom effect. A schematic representation of the ZGB model is the following, 

CO(g)-KS ^ CO(a) 
O2+2S ^ 20(a) 
CO(a) + 0(a) ^ C02(g)+2S 

Here 5* represents an empty site on the surface, g means gas phase and a means adsorbed. For our simulations we 
assume periodic boundary conditions. The time unit is one Monte Carlo Step per Site, MCSS, in which each site is 
visited once, on average. Averages are taken over 3 x 10"* MCSS for each set of parameters after 2 x 10"' warming 
steps. 

III. RESULTS 

Starting from an empty lattice of size lOOx 100, we wait until the system reaches a steady state at a constant partial 
pressure y. We calculate the fraction of sites occupied by CO molecules, the CO coverage (0co) and the O coverage 
{Oq). RcOi is the rate of production of CO2. In Fig. we show the CO coverage, in Fig. [SJb) the O coverage, 
and in Fig. |4] the production rate, for the models analyzed. The first thing to notice is that for the ZGB-n, Hal, and 
Ha2 models, the continuous transition at yi that presents the ZGB model is eliminated: the CO2 production starts 
as soon as y departs from zero (see inset of Fig. 0]). The fact that for y > the surface can not be completely filled 
with O is obvious for the Hal model: the entrance of oxygens requires the existence of empty sites, empty sites that 
can be filled with a CO and then react and create more empty sites. For the Ha2 and ZGB-n models the reason is 
more subtle: each CO-0 reaction produces two empty nn sites, that contrary to what happens in the ZGB model, 
can not be filled by oxygens. The change in the entrance mechanisms for the oxygens also affects the location of the 
discontinuous transition at 7/2 that now has been shifted toward lower values of y. This is an expected result, since it 
is now easier for the surface to be filled by CO. It is not surprising that this effect is particularly strong for the Hal 
model. Notice that, comparing with the ZGB model, for the ZGB-n and the Ha2 models there is a slight increase 
of the maximum value of the CO2 production that occurs very close to the 2/2 transition, while for the Hal model 
the production increases very fast as soon as y departs from zero but the maximum value reached is lower than the 
others. The reactive windows of the ZGB-n and the Ha2 models are larger than the one for the ZGB model. 

Next we include the effect of CO dcsorption by adding a reaction of the form, 

CO(a) ^ CO(g) + S 

Now we perform a CO or O2 adsorption with probability l-fc and a CO desorption with probability k. As we 
already mentioned, this reaction mimics the effect of the temperature in real systems. In Fig. [5] we present the CO 
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coverages for k = 0.02 (a) and k = 0.06 (b). In Fig.|6]we present the O coverages for k = 0.02 (a) and k = 0.06 (b). In 
Fig. [7] we present the CO2 production rates for k — 0.02 (a) and k = 0.06 (b). One effect of including the desorption 
term is that the 1/2 transition is shghtly shifted toward higher values of y. This is expected since the CO-desorption 
reaction makes it more difficult for the surface to be poisoned with CO. Also as can be seen in Fig. IH^a) and (b), 
the production rate increases with the value of k. Notice that the transition at 1/2 seems to become smoother as k 
increases, this is consistent with experimental results and previous simulations that indicate that, for the ZGB model 
with CO desorption, there is a limiting value of k at which the first-order transition line terminates [H. [l9l [20j. Our 
results suggest that, if there is such a critical k for the new models, it is higher than the one for the ZGB. 

IV. CONCLUSIONS 

In this article we have investigated by kinetic Monte Carlo simulations the catalytic oxidation of CO on a surface. 
We have modified the entrance mechanism of the oxygen atoms in the original ZGB model, to take into account 
experimental results that indicate that the oxygens tend to come apart after the adsorption. We explore three 
different entrance mechanisms. In the so-called ZGB-n model the oxygens enter in nnn sites in the surface (see 
Fig. [ij . In the Hal and Ha2 models they enter in the neighboring sites of a selected site that must be empty for the 
Hal and can be filled for the Ha2 (see Fig. [2]). The ZGB model predicts two phase transitions: a continuous one at 
low values of the CO pressure, and a discontinuous one at higher values of y. The continuous transition has not been 
observed experimentally. We found that for all the proposed mechanisms the nonphysical transition disappears and 
the continuous transition is slightly shifted toward lower values of y, the last effect is particularly evident for the Hal 
model. Finally, in order to make the models more realistic, we include a CO desorption term that reproduces the 
effect of the temperature. We found that in all the cases the discontinuous transition is shifted toward higher values 
of y, and that the production rate increases with k. We think that it is possible that the adsorption mechanisms 
presented here occur in different real catalytic processes, and we hope that this study can be useful for identifying 
them. 
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FIG. 1: O2 adsorption mechanisms. One O is adsorbed at the initially empty site i. a) ZGB model: the second O is adsorbed at a nn of 

i. b) ZGB-n model: the second O is adsorbed at a nnn of i. 
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FIG. 2: O2 adsorption mechanisms for the Hal and Ha2 models. The two O atoms arc adsorbed in the vertical or horizontal pair nn to 
i. In the Hal model i must be empty for the adsorption to occur. In the Ha2 model i can be empty, filled with an O, or filled with a CO. 
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FIG. 3: Average values of the CO eoverage (a), and the O coverage (b), shown as functions of the stationary appHed CO pressure y, 

calculated for the different models. 




FIG. 4: CO2 production rate in terms of y for the different models. Notice that for the standard ZGB model the production rate is zero 
until y reaches a minimum value yi !^ 0.38, while for the other models the production starts as soon as iy > as can be seen in the inset. 




FIG. 5: CO coverage when there is a CO desorption probability k. (a) k = 0.02 (b) k = 0.06. 
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FIG. 6: O coverage when there is a CO desorption probability k. (a) k = 0.02 (b) k = 0.06. 
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